Probing the electron-electron interaction in a diffusive gold wire using a controllable 

Josephson junction 
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We have studied the critical current of a diffusive superconductor - normal metal - superconductor 
(SNS) Josephson junction as a function of the electron energy distribution in the normal region. 
This was realized in a 4 terminal device, in which a mesoscopic gold wire between two electron 
reservoirs is coupled in its center to two superconducting electrodes. By varying the length of the 
wire and applying a voltage over it we are able to control the electron distribution function in the 
center of the wire, which forms the normal region of the SNS junction. The observed voltage and 
temperature dependence are in good agreement with the existing theory on diffusive SNS junctions, 
except for low energies. However, an electron-electron interaction time to =10 ps was found, which 
is three orders of magnitude faster than expected from theory. 

PACS numbers: 73.23+r 85.25.Cp 74.50+r 85.25.Am 



The study of electron-electron interaction in metals has 
attracted considerable attention recently. One of the rea- 
sons is an experiment performed by Pothier et al. Q in 
which the electron distribution function in a mesoscopic 
wire connected to large electron reservoirs was measured 
using a superconducting tunnel junction. The shape of 
the distribution function depends, at low temperatures, 
on the effective electron-electron interaction experienced 
by the electrons, and the voltage V c applied over the 
reservoirs. If the wire is short, the electrons will keep 
their energy while traversing the wire, and the distribu- 
tion function will show a double step structure, with a 
separation between the steps of eV c . If the wire is much 
longer, so that the diffusion time rrj through the wire 
strongly exceeds the electron-electron interaction time To, 
then a local thermal equilibrium of the electron system 
will be regained with an effective temperature T e ff, de- 
pending on eV c . The electron interaction time constant 
To obtained from these experiments yields, assuming a lo- 
cal two particle interaction, To ~ 1 ns. for copper, which 
is two orders of magnitudes faster than expected from 
theory §f . 

On the other hand, it is in principle also possible to 
use a diffusive superconductor - normal metal - super- 
conductor junction as a probe to study electron-electron 
interaction effects in a mesoscopic wire. The prediction 
is that J c i? n -product (critical current times normal state 
resistance) of such a SNS junction is very sensitive to the 
exact electron distribution function In this Letter 

we report experiments in which we studied the critical 
current of a SNS junction as a function of the electron 
distribution. This was realized using the device shown 
in Fig. [Ij A diffusive gold wire (control channel) be- 
tween two very large electron reservoirs is coupled to two 



niobium electrodes by means of a cross shaped extension 
in the center of the wire. We study two different device 
geometries, which differ only in the length of the control 
channel. We made 3 devices with a short control chan- 
nel, L contro i = 1 Mm, and two devices with a long control 
channel, L contro i = 9 /im. The behaviour of identical de- 
vices was similar, we therefore present only two devices in 
the remainder of the text: Device 1, with a short control 
channel and device 2 with a long control channel. 




FIG. 1. One of the SNS devices. In this case the length of 
the control channel is 1 fim, The gold control channel (100 nm 
wide) shows a cross in its center, of which the two arms are 
covered by niobium electrodes, thus forming the SNS junction 
in the center of the channel. Contacts at the channel ends 
measure the control voltage. 

The devices are made by means of e-beam lithogra- 
phy using a double layer of PMMA resist and a subse- 
quent lift-off process on top of a thermally oxidised Si 
wafer. The 100 nm wide control channel is evaporated 
first, which consists of 40 nm of gold on top of 5 nm 



1 



Ti, which was used to improve adhesion. Subsequently 
70 nm of niobium was sputter deposited after an in-situ 
argon etching of the gold contacts to ensure a high inter- 
face transparency. In a last step the very thick (475 nm) 
gold reservoirs were deposited. The size of the reservoirs 
is in the order of a millimeter because they should also 
act as effective cooling fins to prevent unwanted electron 
heating pj] , p^ ]. The gold has a diffusion coefficient of 
0.020 m 2 /s which results in an estimated diffusion time 
through the channel of to = 50 ps. and 4 ns. for device 

1 (Lcontroi = IfJ-m) and 2 {L contro i = 9^m) respectively. 
The SNS junctions all have a normal state resistance of 
2.1 f2, and a separation of the niobium electrodes of 375 
nm. 

In the experiment we measured the current-voltage (I-V) 
characteristics of the SNS junction as a function of both 
the voltage V c applied over the control channel as well 
as the bath temperature Tb- In the first case the bath 
temperature is kept at 100 mK to obtain a very sharp 
electron distribution function in the reservoirs. RC fil- 
tering at room temperature and copper powder filtering 
at the bath temperature are used in all measurements to 
reduce external noise and hence, unintentional heating 
of the electrons. From these measurements we obtained 
the critical current. In Fig || we show the results. The 
voltage dependence of device 1, as shown in the bottom 
panel of the figure, shows the transition to a 7r-junction 
at control voltages V c > 0.5 meV, similar to previous ex- 
periments However, the maximum supercurrent in 
the 7r-state is much smaller than expected when assum- 
ing a perfect step distribution function. The top panel of 
Fig. U shows the results obtained in the limit of a ther- 
mal distribution function: The filled circles and squares 
represent the temperature dependence of the I c R n prod- 
uct of device 1 and 2 respectively. The empty squares 
represent the behaviour of the critical current of device 

2 as a function of V c , expressed in terms of the effective 
temperature T e ff, assuming a perfect thermalisation of 
the electrons in the wire . The effective electron tem- 
perature in the center of the control channel can now be 
calculated using the Wiedemann- Franz law M: 



T, 



eff 



= \ T B 2 + (a-Vy 



(1) 



with a=3.2 K/mV. All curves show the expected mono- 
tonic decrease in I c with increased temperature ]To[ | , and 
at higher temperatures all curves lay essentially on top 
of each other, with the exception that I c R n vs. T e ff 
(open squares) is somewhat larger than I c R n vs. Tb 
(solid squares) at higher temperatures. This indicates 
that T e ff increases slower with increasing V c than ex- 
pected, which is probably due to the onset of electron - 
phonon interactions at higher electron energies [jl4| . This 
leads to an extra cooling of the electrons, and thus to a 
larger supercurrent. A very striking difference between 
both devices, which is counterintuitive and seems to be 



in disagreement with theoretical predictions JsJ , arises at 
low temperatures, T < 800 mK.: device 1, with the lowest 
control channel length, and thus the strongest coupling to 
the normal reservoirs, shows a much higher 7 c i?„-product 
(87 /ieV), than device 2 (I c R n =56 /xeV), whereas one 
would expect a reduction in supercurrent with increased 
coupling to a normal reservoir. This observation is con- 
sistent with a previous experiment, which yielded an even 
higher 7 c i? n -product for a system with an even stronger 
coupling to the normal reservoirs |pf . 
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FIG. 2. I c R n as a function of temperature and voltage for 
both devices. The top panel shows the temperature depen- 
dence of device 1 (circles) and device 2 (filled squares), and 
the voltage dependence of device 2 (open squares) expressed 
in units of T e / / . The bottom shows the voltage dependence 
of device 1. The best fits to the theory are given by the solid 
lines, assuming a distribution function as shown in the insets. 

The next step in the analysis is to compare the data 
to the existing theory on diffusive SNS junctions. This 
theory predicts that supercurrent will be carried by a su- 
percurrent carrying density of states, Im(J(e)), which can 
be calculated using the quasi-classical Green's function 
theory The positive and negative parts of the su- 

percurrent carrying density of states, as shown in of Fig. 
||, represent, at a given phase difference 7r/2, energy de- 
pendent contributions to the supercurrent in the positive 
and negative direction. The critical supercurrent now 
depends strongly on the occupation of this continuum of 
states, and thus on the electron distribution function in 
the normal region, and can be calculated according to 
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I C R N = / de[l - 2/(e)]/m(J(e)) (2) 

J — oo 

where ImJ((e) is the supercurrent carrying density of 
states and f(e) the electron distribution function. The 
transition to a 7r-junction in case of a perfect step func- 
tion can be understood from this, for in this case the 
distribution function will be exactly 0.5 over a region eV c 
around the Fermi energy, resulting in a zero contribu- 
tion to the supercurrent over this energy range. At a 
large enough value of V c , its magnitude depending on 
Eth-, all positive contributions to the supercurrent will 
be blocked, which obviously changes the direction of the 
supercurrent. 




FIG. 3. Supercurrent carrying density of states, using 
A/E th =70, A=1.5 meV (Niobium) and E th =2l^,eV. (from 
ref. [7,8]) 

We tried to fit the temperature and voltage dependence 
of the J c i? n -product, using eqn. ^| in which the Thou- 
less energy as well as the exact shape of the distribution 
function were used as fit parameters. However, the shape 
of the distribution function can be calculated as a func- 
tion of the effective electron-electron interaction time in 
the control channel, which leaves only E t h and To as fit 
parameters. 

The calculation of the distribution function was per- 
formed using the method described in Q, assuming a 
local two particle interaction with an interaction kernel 
given by -^pz, with e the transferred energy and t the 
electron-electron interaction time constant. In this model 
the shape of the distribution function depends on 2 pa- 
rameters: The ratio of the electron interaction time with 
the electron diffusion time, tq/tb and on kT e /eV Cl which 
is the ratio of the electron temperature in the reservoirs 
T e and the control voltage applied. T e in the reservoirs 
cannot be taken as a constant and equal to Tb, because 
the electron heat conductance of the reservoirs is finite 
and electron-phonon interaction is strongly reduced at 
low temperatures. These arguments combined with the 
fact that rather large voltages are used in these experi- 



ments, make it necessary that the heating of the electrons 
in the reservoirs is calculated explicitly. This was done 
using a model presented by Henny et al. Jll],[l2|]. The 
result is that, at Tb ~ 100 mK, T e increases with in- 
creasing control voltages, but this increase is linear if V c 
> 0.4 mV. In this limit kT e /eV c is constant, which im- 
plies that the shape of the distribution function is also 
constant. Higher values of the applied voltage merely 
change the energy range of the distribution function, not 
its shape 0]. 

The results of the fits are shown by the solid lines in Fig. 

in which the maximum value of the J c i?„-product of 
device 1 was used as normalisation constant. The inset 
shows the distribution functions (in units of E/eV c ) in 
the limit where kT/eV c of the reservoirs is constant. We 
found a value of the Thouless energy given by Eth—Tl 
meV, which is in good agreement with the expected I c R n - 
product of a diffusive SNS junction, = n [pli|| , as 

well as the value reported in || for the cross geometry 
discussed here, I f ? Rn = 5. The effective Nb electrode 
separation using this value of E t h is l=800nm, which is 
in between the minimum Nb separation and the maxi- 
mum extent of the gold under the Nb. It is clear from 
Fig ||that the agreement between the experiment and the 
calculations is excellent, apart from the low voltage re- 
gion of device 2 discussed previously. However, from the 
insets is is also clear that the shape of the distribution 
functions indicates a significant electron-electron interac- 
tion: The distribution function has a Fermi-Dirac shape 
in the case of device 2, which was expected, but it is also 
very rounded already in the case of device 1 . This strong 
rounding is responsible for the small magnitude of the su- 
percurrent in the 7r state. The interaction time constant 
ro obtained from this analysis is device independent, and 
given by To ~ 10 ps. We have found that choosing a 
larger relaxation time constant in combination with an- 
other normalisation constant or another value of E t h does 
not yield reasonable and consistent fits for both devices. 
Moreover, the dependence of the 7 c i?„-product on the 
control voltage is extremely sensitive to the effective elec- 
tron relaxation in the control channel, tb/tq, and thus 
on the exact shape of the distribution function, as shown 
at the bottom of Fig. ||. Here the calculated distribution 
function as well as the resulting I c R n - V c behaviour is 
plotted for 5 different distribution functions, using five 
different ratio's of tq/tb- It is clear that even small de- 
viations from a thermal distribution function (curve b) 
results in a quite strong difference in the voltage depen- 
dence of the I c i?„-product, even leading to the transition 
to a 7r-junction at relatively rounded distribution func- 
tions (curve c). 

Up till now we have based our analysis on the assump- 
tion that Im(J(e)) is correct. However, the small values 
of the 7r supercurrent could in principle also be explained 
assuming a perfect step distribution function if the nega- 
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tive contributions to the supercurrent of Im(J(e)) would 
be smaller. This assumption yields the fit of the I c R n 
vs. V c data of device 1 as shown in panel A of Fig. 
U However, if this form of 7m(J(e)) is used to calculate 
the temperature behaviour, we get a strong disagreement 
between measurements and theory, as shown in panel B. 
We therefore conclude that the theory which describes 
the supercurrent as a function of the distribution func- 
tion is in principle correct, except at low energies, for it 
fails to predict the difference in I c i? ra -product between 
the devices at low temperatures. 
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FIG. 4. The two top figures show a theoretical fit to 
both the voltage and temperature dependence of device 2, 
assuming a perfect step distribution function and, as a result, 
smaller negative contributions in Im(J(e)). The bottom figure 
shows the strong dependence of the I c R n vs. V c behaviour on 
the effective electron-electron interaction in the control chan- 
nel, to/td- The inset shows the corresponding distribution 
functions. t /t d = : a: 0, b 0.1 c 0.2 d 0.4 e 1. 

On the other hand, the value of To implies an upper 
bound for the quasiparticle lifetime, valid over the energy 
range of these experiments (0.01 - 1.8 meV), given by 
r(E) < T /ln(E/0.01 meV) §. This yields a maximum 
quasiparticle lifetime, at E=0.01 meV, given by T e _ e ~ 
10 ps. However, the electron out sacttering time can be 
calculated using the existing theory on electron-electron 
interactions in a diffusive ID wire, according to 



S the area cross section of the wire. However, this yields, 
using n =1.9-10 28 m^eV' 1 r ee (0.01 meV)=80 ns, a dif- 
ference of three orders of magnitude with our result. We 
note that experiments on simular wires using supercon- 
ducting tunnel junctions [jl 16 give a a value of tq 1 ns 
and 0.1 ns for copper and gold respectively. The differ- 
ence between our data and these measurements can be 
explained by the fact that we used a thin Ti adhesion 
layer (which has a very fast phase relaxation at low tem- 
peratures, caused by electron-electron interaction |l9[l ). 
Alternatively the strong coupling between the niobium 
and the gold in the junction might influence the electron 
interaction. However, the difference with the existing 
theory is more than 3 orders of magnitude, which 
can hardly be explained using such arguments. Very re- 
cently JlTj the experimental procedure of Ref. Q| to ob- 
tain the electron distribution was questioned, as well as 
the method discussed here The prediction is that 

the shot noise in the control channel associated with the 
non-equilibrium distribution function yields a smearing 
of the density of states in the superconducting electrodes, 
resulting in a supercurrent reduction. Our observations 
are not in agreement with the predicted Jl8| behaviour 
of I c vs. V c (see bottom panel of Fig. ||). Also, assuming 
the theoretical value To ss 80 ns would result in a tran- 
sition to a 7r-state of device 2 as well, which we did not 
observe in the experiment. 
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(E) = V2Hn S 




(3) 



with n the density of states around the Fermi level, and 
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